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ABSTRACT 

The three-dimensional (3D) structure of solar active region NOAA 7986 observed on 
1996 August 30 with the Extrem-ultraviolet Imaging Telescope (EIT) onboard the Solar 
and Heliospheric Observatory (SoHO) is analyzed. We develop a new method of Dynamic 
Stereoscopy to reconstruct the 3D geometry of dynamically changing loops, which allows us 
to determine the orientation of the loop plane with respect to the line-of-sight, a prerequisite 
to correct properly for projection effects in 3D loop models. With this method and the 
filter-ratio technique applied to EIT 171 A and 195 A images we determine the 3D coordinates 
[x(s), t/(s), z(s)], the loop width w(s), the electron density n e (s), and the electron temperature 
T e (s) as function of the loop length s for 30 loop segments. Fitting the loop densities with 
an exponential density model n e (A) we find that the so inferred scale height temperatures, 

T e A = 1.22 ± 0.23 MK, match closely the EIT filter-ratio temperatures, Tf IT — 1.21 ± 0.06 MK. 

We conclude that these rather large-scale loops (with heights of h w 50 — 200 Mm) that dominate 
EIT 171 A images are close to thermal equilibrium. Most of the loops show no significant 
thickness variation u?(s), but many exhibit a trend of increasing temperature (dT/ds > 0) above 
the foot point. 


Subject headings: Sun: corona — Sun : X-rays, gamma rays — Sun: EUV 


1. INTRODUCTION 
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2. STEREOSCOPIC DATA ANALYSIS 
2.1. Data Set 

The investigated active region is a long-lived coronal structure that was present during several solar 
rotations, from its apparition in July 1996 until its disappearance in September 1996 (Hudson et, al. 1998; 
Harvey & Hudson 1998; Portier-Fozzani et al.1998), numbered as NOAA 7978, 7981, 7986 during subsequent 
rotations. We concentrate here on the central meridian transit on 1996 Aug 30, when the dipolar magnetic 
field structure offered the most favorable perspective to disentangle the “jungle” of nested loops. 

An Fe IX/X image recorded with SOHO/EIT (Delaboudiniere et al. 1995) at a wavelength of 171 A on 
1996-Aug-30, 0020:14 UT, is shown in Fig.l (top). For stereoscopic correlations we will use also EIT images 
from the previous ( 1996-Aug-29, 0015:15 UT) and following day (1996-Aug-31, 0010:14 UT). The multi-loop 
structure of this active region is clearer visible in the high-pass filtered rendering shown in Fig.l (bottom) 
The filtered image is simply created by subtracting a smoothed image (with a boxcar average over 3x3 
pixels) from the original image. The original (full disk) image has a pixel size of 2.585 and was recorded 
with an exposure time of 3.5 s. The absolute coordinate system of the full-disk image was established by 
fitting a circle to the solar limb (at 30 limb positions). The accuracy of the so-defined Sun center position 
is estimated to be (T % 1 pixel/v/30 « 0.2 pixel. The offset of the Sun center position (x' 0 ,?/o) provided 
(by an automatic limb fitting routine) in the FITS header of the archive EIT image with respect to our 
value (a?o> Vo) is found to be (x f 0 — £o) = +0.5 and (yj — yo) = —3.9 pixels. For the solar radius we find a 
difference of (r{, - r 0 ) = —1.2 pixels, yielding a pixel size and uncertainty of Ax = 2.585” ± 0.004”. 


2.2. Dynamic Stereoscopy Method 

In order to analyze the 3-dimensional structure of coronal loops we develop a new technique we might 
call “Dynamic Stereoscopy” } opposed to u Static Stereoscopy ” , where the solar rotation is used to vary the 
aspect angle of otherwise static structures (e.g. Loughhead, Chen & Wang 1984; Berton fe Sakurai 1985; 
Aschwanden & Bastian 1994a, b; Aschwanden et al. 1995). The innovative feature of this new technique is 
that spatial structures, e.g. coronal loops, are allowed to evolve dynamically during the time interval over 
which the stereoscopic correlation is performed. 

In the Dynamic Stereoscopy method we take advantage of the fact that, the global magnetic field 
is slowly evolving (say during a day) compared with heating and cooling processes in coronal loops. 
Consequently, the coronal magnetic field B(x,t) can be considered as invariant over short time scales, while 
the confined plasma can flow through “magnetic conduits” in a highly dynamic manner. If a specific coronal 
fluxtube following a field line B(x,U) is loaded with bright plasma at time U, the same fluxtube may be 
cooled down at time to (say a day later) and invisible at the same observed wavelength, while heating may 
occur in an adjacent fluxtube B(x + Ax,h), which was dark at time and appears now bright at time to- 
For adjacent fluxtubes, the two field lines B(x) and B(x + Ax) will run almost parallel, a property we will 
exploit in our dynamic stereoscopy method. Our method is applicable to coronal structures that meet the 
following two conditions: 

T The global magnetic vector field B(x } t) is static (or slowly-varying) during the time interval over 
which stereoscopic correlations are performed (typically I day). The magnetic field can be traced out 
by optically-thin emission (e.g. in soft X-ray or EUV wavelengths). 
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2. At least one footpoint of an observed coronal loop is identifiable, which can be used as a reference 
level of the altitude. For EUV emission we assume that the altitude of a loop footpoint is located in 
the lower corona above the chromosphere, at an altitude of hj 00t & 2500 km above the photosphere. 


We outline briefly the numerical procedure of our implementation of the dynamic stereoscopy method, 
while the mathematical coordinate transformations are given in Appendix A. The projected geometry of 
a loop segment in an image at time t i is traced out by a series of image coordinates (£*,«/*)>* = 1 
starting at footpoint position (xi,yi), assumed to be anchored at height hf oot (Fig. 2). Two additional 
variables to characterize the 3D geometry of the loop segment are the azimuth angle a of the footpoint 
baseline, and the inclination angle d of the mean loop plane (intersecting the footpoint baseline, see Fig. 2). 
The procedure of stereoscopic correlation (illustrated in Fig. 3) includes the following steps: 


1 . Measuring of positions (x iy yi),i = 1 , by tracing a loop segment in an image recorded at time 1 1 , 

starting with the primary footpoint at (ar i , i/i) ■ 

2. Estimating the position of the secondary footpoint (x fo, Vf?) to obtain the azimuth angle a of the 

footpoint baseline. If the full length of the loop can be traced, the secondary footpoint is just given by 
the last point x n ,y ny and the tangens of the azimuth angle a corresponds to the ratio of the latitude 
(bp 2 — ) and longitude difference (/ f 2 — /fJ of the two footpoints, i.e. 


tan(a) 


(M - b Fl ) 

(If 2 - hi) 


( 1 ) 


However, most of the loops analyzed here extend over several density scale heights, so that only a 
segment above the primary footpoint can be traced out, while the location of the secondary footpoint 
has to be estimated based on the general dipole characteristic of the ambient magnetic topology. 

We repeat the localization of the secondary footpoint 5 times and obtain from these estimates 
aj , j 1 , 5 a mean and standard deviation a db da. 

3. The loop positions (jr^y,-) measured in the image at time are then transformed into heliographic 
longitude and latitude coordinates (/,j , bij) and altitudes (/i 2 j)i based on the azimuth angle a 

of the footpoint baseline and the variable inclination angle i9j, which is varied over a range of 
—90° < flj < +90° in increments of = 1°. 

4. The heliographic coordinate is then transformed to the time t 2 of the second stereoscopic pair 

image, /{•(* 2 ), using the solar differential rotation rate applied to the time interval (t 2 - fi). We use 
the differential rotation rate specified by Allen (1973), 


l'(t 2 ) = 1(h) + [13.45° + 3°sin 2 (B)] ^ . (2) 

The heliographic latitude bij and altitude hij are assumed to be constant during the considered time 
interval. 


5. In the stereoscopic pair image at time t 2 we calculate the image coordinates {x 1 ^, y' 2j ) of the projected 
loop structure. Parallel to these loop curves (with typical lengths of n, = 50 — 200 pixels) we extract 
image stripes of some width (n w = 16 pixels) by interpolating the image brightness F(x, y) at the 
positions of the curved coordinate grid. 
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6. The stretched 2D image stripes (n s x n w pixels) are then scanned for parallel brightness ridges, caused 
by “dynamic” structures that are co-aligned (or parallel-displaced) to the loop projection in image t*i 
(for illustration see examples shown in lower panels of Fig. 3). This scanning process is numerically 
performed by measuring the total length L(d^) of coaligned brightness ridges that are detected in 
each stripe k associated with each angle dj. The maximum of this quantity (rnax[L(d k)]) 18 used to 
infer the most likely value of the inclination angle d of the loop plane. 

7. The same procedure is repeated in forward and backward direction in time. The mean and standard 
deviation of d ± dd is determined by averaging the two stereoscopic solutions (d=L day). 

The independent stereoscopic correlation in forward (-f 1 day) and backward ( — 1 day) direction provides a 
useful redundancy of the solution. The time difference of ±1 day corresponds to an aspect angle change of 
±13.5°. Except for steps 1-2, which constitute the definition of a selected loop feature, all other steps (3-7) 
of the stereoscopic correlation are performed automatically by a numeric code without human interaction. 
The determination of the loop inclination angle d is therefore achieved in a most objective way, within the 
principle of dynamic stereoscopy. 


2.3. Loop Geometry 

With the dynamic stereoscopy method described above we analyzed the 3-dimensional (3D) coordinates 
of 30 loops from the EIT 171 A image on 1996-Aug-30 (Fig.4, middle column). The true 3D geometry 
of (the central axis) of a coronal loop can be characterized with three orthogonal spatial coordinates 
[x(si),y(si), z(si)],i = 1, n, parametrized by the loop length parameter s*. When we trace a loop 
structure in an image (see Fig.4, panels in middle column), we can accurately measure the two coordinates 
[x(sj), y(si )], without imposing any geometric constraint on its shape, opposed to a method that fits a 
predefined geometric model (e.g. an elliptic projection of a circular geometry). We impose only some 
constraints on the third coordinate z(s,), namely assuming that the loop segment is confined into a plane, 
whose orientation we quantify with two free parameters (Fig. 2): with the azimuth angle a (of the loop 
footpoint baseline) and the inclination angle d (with respect to the vertical). One additional constraint is 
also introduced by the reference level h(s\) of the first footpoint position, assumed to be located at a fixed 
height of h($\) = hf 00 t ^ 2.5 Mm. 

Some geometric elements of the analyzed 30 loops are listed in Table 1: the heliographic longitude 
(/i) and latitude (bi) position of the primary footpoint [x(si), y{s \ ), h{s \ ) — hf oot ], the azimuth angle a of 
the footpoint baseline measured at the primary footpoint, and the inclination angle d of the loop plane. 
The average heliographic position of the 30 loop footpoints is < lp x >= 251.0°, < bf' l >= —11.8°, which is 
slightly southward of the Sun center position at this time, lo = 255.8°, bo — 7.2°. 

The average azimuth angle (modulo 180°) of the 30 loop footpoint baselines is a ~ —3° ± 10°, which 
represents the global orientation of the large-scale dipolar magnetic field that dominates the active region, 
that was used as a guide to estimate the azimuth angle of the footpoint baseline for indiviudal loop 
segments. The only complete loop that could be traced out (without gaps between the footpoints) is Loop 
#1, which has an azimuth angle of a = 15° ± 1°. 

The inclination angles i? of the loop planes cover a large range from d = -56° to d — 4-69°, having 
an average of < d >= 7° ± 37°. The southern loops (Loops #1-11, 23-30) show all an inclination toward 
south (with d negative if the primary footpoint is to the East), while the northern loops (Loops #12-22) 
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show a systematic inclination toward North (with d positive if the primary footpoint is to the East). This 
fan-shaped divergence of loop planes is consistent with the overall magnetic dipolar field, having the dipole 
axis aligned to the East- West direction. 

We visualize the 3D structure of the stereoscopically reconstructed 30 loops in Fig. 5, where different 
perspectives and viewing angles are displayed. The traced segments (Fig. 4) of the reconstructed loops 
are marked with thick lines in Fig. 5, while the thin lines represent circular segments in the loop plane, 
constrained by the two footpoints and the endpoint of the traced segment. We rotated these reconstructed 
30 loops by -7.2 days to the East (Fig. 5 bottom left), in order to illustrate the distribution of inclination 
angles. The group of loops that are inclined to the South in our EIT image of 1996- Aug-30, are also found 
to have a similar configuration (with similar loop heights and loop plane inclinations) in an EIT image 
observed 7.2 days earlier, when this active region crosses the East limb (Fig. 5 top left). 


2.4. Loop Background Subtraction 

We parametrize the positions of the traced loop segments from the image coordinates [a: («,•), y{s , )] 
as function of the (projected) loop length parameter Sj, by interpolating the coordinates with a constant 
resolution of As = s t+ i — s a = 1 pixel in the image plane. These positions mark the central axes of 
the analyzed loops. For single-loop analysis it is convenient to introduce a coordinate grid [s{ y tj] that 
is stretched out along the loop axis, with tj a coordinate orthogonal to the loop axis. The projections 
[r(S| , tj), y(si, tj)] of these curved coordinate grids are shown in Fig. 6 (top right panel). We parameterized 
both coordinates [s i} tj] with a uniform resolution of 1 pixel, and have chosen a width of 16 pixels for the 
width of the stripes (2j), symmetrically bracketing the central loop axis. We show the flux F(tj) as function 
of the loop cross-section tj for each loop (#1—30), and for each position s* along the central loop axis with 
a step of Asi = 1 pixel in Fig. 6, measured from the EIT 171 A image of 1996-Aug-30. 

In a next step we attempt to separate the loop- associated fluxes from the loop-unrelated background. 
This is a very crucial step to determine the correct emission measure and electron density in a given loop. 
This task is difficult because most of the loops are very closely spaced and separated only by a few pixels 
at their primary footpoint (see Fig. 6). Very few loops occur in an isolated environment (e.g. such as 
loop #25, see Figs.4-6). For many cross-sections there is not enough separation between adjacent loops 
to model the loop- unrelated background properly. The fact that the flux-unrelated background makes 
up typically 50-90% of the total EUV flux measured at the same position (see Fig.7), indicates that we 
can separate out only a fraction of super-imposed and nested loops, like the top-most elements on the 
topological surface of a “strang of spaghettis”. We tested various methods and found the following to 
be least susceptible to confusion by adjacent loops. We calculated the background profile F B (tj) to the 
observed flux F(tj) by cubic spline interpolation between various intervals [<i,<2]> which were varied over a 
range of t 0 - t\ = 1...4 and t 2 - to = 1...4 pixels (or 2-7 Mm for the half width of the loop cross-section) 
on both sides of the central loop axis (at * 0 ). From the varied loop boundaries [t u t 2 ] we select those values 
for the best background estimate that maximize the flux of the loop cross-section, i.e. which yield the 
maximum of /[^(fj) — F B (tj)]dtj. This method has the advantage to adjust for loop thickness variations, 
for offsets in tracing of the central axis, and for co-alignment errors between the 171 and 195 A image in 
the use of the filter-ratio technique. The so-determined loop-associated fluxes are shown with grey areas for 
each cross-section F{tj)\ s=Si in Fig. 6. The results show that the allowed loop half-width range of ±1 - 4 
pixels separates most of the loops reasonably, except for occasional double loop detections near the primary 
footpoint. Such loop segments where the loop separation fails will be excluded in further analysis. 
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2.5. Loop Cross-Sections 

We measured the loop width u;(s) as function of the loop length parameter s, using the definition of 
the equivalent width w(s), 

/ F(s,tj) — F j b(s,<j) dtj /g\ 

U S max[F(s,tj) — Fb(sJj)] 

These loop widths w(s) are shown as function of the loop length s in Fig. 8 for the 10 loops that are least 
confused by adjacent loops, as it can be judged from the cross-sections F(tj) | J=Jl shown in Fig. 6 (Loop 
#1, 8, 11, 14, 15, 19, 20, 21, 25, 28). Performing a linear regression fit to the observed values u>(s) we 
find a significant variation of the loop thickness only for two of them (Loop ^20 and 28). To quantify the 
variation of the loop thickness we calculated a loop divergence factor , defined by the average width in the 
upper part (s max / 2 < s < s max ) to the lower part (0 < s < s max / 2) of the traced loop segment. We remind 
that the traced loop segments generally extend over about one density scale height, but do often not reach 
the loop top (except for the smallest loop, #1). The loop divergence factors and their uncertainties are 
shown in Fig. 8 (bottom right panel) for each loop. We caution that some of the loop thickness variation 
near the footpoints is due to separation problems of closely-spaced adjacent loops (as it can be judged from 
Fig. 6). A histogram of average loop widths is shown in Fig. 8 (top right panel), while the individual values 
u> and their mean and standard deviation (te = 7100 ± 800 km) is also listed in Table 1. The preference 
for such a narrow range of loop diameters is perhaps an instrumental resolution bias, because the finest 
recognizable structures are most likely to be seen at a scale corresponding to the size of a few pixels. 


2.6. Loop Densities and Scale Heights 


For electron density and temperature diagnostics we are using a filter ratio technique applied to the 
FIT 171 A and 195 A wavelength images, based on the most recent EIT standard software (status of 
February 1989). The resulting emission measures EM and temperatures Tf IT are based on the calculation 
of synthetic spectra using the CHIANTI data base, containing some 1400 emission lines in the 150-400 A 
wavelength range. Comparing these current standard values (as referred to in this study) with recent new 
calibration data obtained from a rocket flight with a duplicate EIT instrument (led by Dan Moses, NEL), 
the emission measure values quoted here would have to be corrected by ps 20%, or the electron density by 
& —40%, respectively, while the temperature values remain basically unchanged. 

To determine the electron density n e (s) along individual loops, we use the background-subtracted EIT 
fluxes Fi oop (s) = F{s) - F B (s) in the filter ratios and the loop widths io(s). An additional important loop 
parameter is the line-of-sight angle ^ ( s ) , which provides a correction factor of the effective column depth for 
a loop with circular cross-section specified by a diameter tc(s), i.e. u; 2 (s) = w(s)/ cos[i/’(s)] (see Appendix 
B). With this parameterization we define the density n e (s) along a loop by 


n e (s) 


lEM(s) EM (s) cos[^(s)] 


«h(s) 


tc(s) 


( 4 ) 


with the loop length s(x,y) parametrized as a function of the image position (x,t/), from where the EIT 
emission measure EM(x,y) is measured. Because the line-of-sight angle ip(s) is very sensitive to the loop 
orientation, correct values of the electron density n e (s) can only be obtained from an appropriate 3D model 
of the loop (constrained by stereoscopic correlations here). The projection effect of the loop curvature on 
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the effective column depth w z (s), and the effect of the inclination angle $ of the loop plane on the inferred 
density scale height A(tf) are illustrated in Fig. 9. 

The electron density n e (s) calculated from Eq.4 is shown graphically for the 10 least-confused loops (# 
1,8,11,14,15,19,20,21,25,28) in Fig. 10 (left panels). Because the height dependence s(h) of the loop length is 
known from our stereoscopic reconstruction (displayed in Fig. 5), we can directly obtain the parametrization 
n e [s(/i)] »-»■ n e (h) and fit an exponential density model, 

n e {h) = n e0 exp(- ^ ) (5) 


to obtain a scale-height temperature T e \ which is defined by (e.g. Lang 1980, p.285) 


A(Te) = 


ksT e 

jj. m H g 


46 <ij [Mmi ' 


( 6 ) 


with k B the Boltzmann constant, /< the mean molecular weight (// 1.4 for the solar corona), the 

mass of the hydrogen atom, and g the acceleration of gravity at the solar surface. The so obtained scale 
heights A, with a mean of A = 55 ± 10 Mm, and the inferred scale height temperatures T A , with a mean of 
T A = 1.22 ± 0.23, are listed in Table 1 for each of the analyzed 30 loops. Loop segment ranges that are 
obviously confused by adjacent or crossing loops (as it can be judged from Fig. 6), have been excluded in the 
fitting of the scale height model. We find that most of the analyzed loop segments fit closely an exponential 
density model (see Fig. 10 left). Deviations from an exponential density model can often be explained by 
uncertainties in the background subtraction, or by confusion from adjacent or overlying loops. 


2.7. Loop Temperatures 

Independently of the scale height temperature T A , we can also determine the temperature directly 
from the EIT filter ratio, which moreover provides a temperature differentiation along the loop, Tf lT (s ). 
Since our loop definitions are based on tracing of EIT 171 A image, we use only the filter ratio of EIT 
171 A (Fe IX, X) and 195 A (Fe XII), which is sensitive in the temperature regime of T e = 1.0 - 1.5 
MK. We are using the spatial loop definition [ar(s),y(s)] based on the 171 A image, and apply the same 
background-subtraction technique to the 195 A image, before we determine the temperature from the filter 
ratio (F 171 - F^ 71 )/(F 195 - F^ 95 ). Because our background-subtraction technique has some tolerance (of 
±1-4 pixels) in the localization of the loop cross-section (by maximizing the flux integrated over the loop 
cross-section), the filter-ratio is not susceptible to small co- alignment errors between the 171 and 195 A 
image. The employed background-subtraction technique requires also a correlated structure (with a width 
of 2-8 pixels) in both wavelengths, while larger or diffuse structures with possibly different temperatures are 
safely subtracted out. 

The filter-ratio temperatures Tf lT averaged over the loop segments are tabulated in Table 1, with a 
mean of Tf IT = 1.21 ±0.06 MK. The distribution of filter-ratio temperatures N(Tf IT ) is shown in Fig. 11 
(bottom left), along with the distribution of scale height temperatures N ( F A ) (Fig. 11 top left), which have 
almost identical means. The range of scale height temperatures (±0.23 MK) is broader than the range EIT 
filter-ratio temperatures (±0.06 MK), probably because of systematic errors in background subtraction and 
loop separation. This is also consistent with the scatterplot of the two temperature definitions (Fig. 11 top 
right), where no obvious correlation is seen. Despite of these unavoidable uncertainties in the background 
subtraction, it is remarkable that the means of the two independently determined temperatures conincides 
so closely. 
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2.8. Loop Temperature Gradients 

The detailed variation of the temperature Tf fT (s) along the loop length s is shown for the 10 
least-confused loops in Fig. 10 (right panels). We note that the filter-ratio temperature varies sometimes 
discontinously along the loop, e.g. there is a jump from T e = 1.35 MK to T e — 1.1 MK at s = 70 Mm in 
Loop #1 (Fig. 10, top right), which may be caused by contamination from a hotter loop that is located 
almost parallel to Loop #1 at s < 70 Mm (see cross-sections in Fig. 6). Such confusion problems can only 
be identified in hindsight. Despite of such confusion problems, there seems to be a trend of a positive 
temperature gradient dT/ds > 0 above the footpoint. To estimate these average temperature gradients 
(without correcting for multi-loop confusion) we performed a linear regression T e (s) for all loops. Most 
significant temperature gradients are found for Loop # 11 (dT/ds = +0.0049 MK/Mm), for Loop #20 
(dT/ds = +0.0012 MK/Mm), and for Loop #25 (dT/ds = +0.0054 MK/Mm). The distribution of 
temperature gradients N (dT/ds) is shown in Fig. 11 (bottom right), revealing that % 75% of the loops have 
a positive temperature gradient dT/ds > 0 within the first scale height above their foot points. Higher parts 
(h > 1A) of these loops are not detectable in EIT images due unsufficient density contrast (one scale height 
corresponds to a factor of w 3 in density or a factor of 10 in emission measure or EIT flux). 


3. PHYSICAL INTERPRETATION AND DISCUSSION 
4. CONCLUSIONS 
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A. (APPENDIX) Heliographic Coordinate Systems and Transformations 

For analysis of observed images, for time-dependent coordinate transformations that take the solar 
rotation into account (as needed in stereoscopic correlations), and for convenient definitions of loop 
geometries we define 3 different coordinate systems: 

Image Coordinate System fa,y,z) : The (x,y) coordinates refer to the x-axis and y-axis of an observed 
image, while the coordinate ( 2 ) is orthogonal to the image, or parallel to the line-of-sight direction, defined 
positive toward the observer. The origin (x,y, z) = (0,0,0) of this coordinate system is most conveniently 
assumed at the Sun center position. A solar FITS image should contain the position of the Sun center in 
pixel units (i x0 J y o) (in FITS header CRPIX1, CRP1X2 or E.XCEN, E.YCEN ), the pixel size (Ax, Ay) 
in units or arcseconds (in FITS header CDELT1, CDELT2) , and the solar radius i r 0 in pixel units (in 
FITS header SOLAR-R, or E.XSMD, E-YSMD if the semi-diameters of an ellipse are fitted). With this 
information, a pixel position (i, j) can then be converted into the coordinate system (x,y) by 

x, = Ax (i — i x o) (Al) 

Vi - Ay (j - .7*0) , ( A2 ) 


where Ax=arcseconds/pixel for (x,y) in units of arcseconds, or Ax = it@/* r o, with Rq — 696,000 kin, if 
physical length units (km) are preferred. 

Heliographic Coordinate System (l,b,r) : The heliographic coordinate system is co-rotating with the 
solar surface. A position on the solar surface is generally specified by heliographic longitude and latitude 
coordiantes (/,&) (in units of heliographic degrees), with reference to the Carrington rotation grid. The 
heliographic longitude and latitude [/o(i),6o(01 Sun cen f er an d the position angle P(t) of the solar 

rotation axis for a given time t are published in The Astronomical Almanach (Nautical Almanac Office, 
NRL, Washington DC). The 2-dimensional spherical coordinate system (/,&) can be generalized into a 
3-dimensional coordinate system by incorporating the height h above the solar surface, which can be 
expressed as a dimensionless distance to the Sun center (in units of solar radii), 


, it V 

r = (! + -5“) ■ 


R, 


'© 


(A3) 


The transformation from the 3D heliographic coordinate system (l, b, r) into image coordinates (x , y, z) can 
be accomplished by applying a series of 4 rotations to the (normalized) vector (0, 0, r), 


/ x/Rq \ 

v/ r q 

\ z/R& / 


/ COS (P + Po) 
j sin(P + Fa) 

V 0 

/ cos(/o — /) 0 

0 1 

\ sin(/o — l) 0 


-sin(F + F 0 ) 0 

cos (F + Fo) 0 

0 1 

— sin(/ 0 - /) \ 

0 

cos(/ 0 - /) / 


/ 1 0 0 

0 cos(6 0 ) — sin(6 0 ) 

\ 0 sin(6 0 ) cos(6 0 ) 


1 0 0 
0 cos(— b) — sin(— 6) 
0 sin(— b) cos(— b) 


0 

0 


r 


(A4) 


where (/oi^o) a re the heliographic longitude and latitude of the Sun center, P is the position angle of the 
solar rotation axis with respect to the North-South direction (defined positive towards East), and P 0 is the 
image rotation (roll) angle with respect to the North-South direction, (P + Pq = 0 for images rotated to 
solar North). In stereoscopic correlations, only the longitude of the Sun center, lo(t), is time-dependent in 
first order (according to the solar rotation rate), while b 0 (t) and P(t) are slowly-varying, and thus almost 
constant for short tirfie intervals. 
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Loop Plane Coordinate System (X,Y,Z) : To parametrize coronal loops it is also convenient to introduce 
a cartesian system that is aligned with the loop footpoint baseline (A -axis) and coincides with the loop plane 
(X— Z plane, Y = 0). The transformation of loop coordinates (X, Y = 0, Z ) into cartesian coordinate system 
(X',Y',Z') that is aligned with the heliographic coordinate system (/, 6, r) can simply be accomplished 
with help of two rotations, 

/ X* \ / cos(o) — sin(a) 0 \ / 1 0 0 \ / A \ 

I Y f I = J sin(a) cos(a) 0 j I 0 cos(t?) sin(i?) 11^1, (A5) 

\ Z f / \ 0 0 1 / \ 0 -sin(tf) cos(i?) } \ Z ) 

where the azimuth angle a denotes the angle between the loop footpoint baseline and the East- West 
direction, and represents the inclination or tilt angle between the loop plane and the vertical to the solar 
surface (Fig. 2). Placing the origin of the loop coordinate system [X — 0, Y = 0, Z = 0] (which is also 
the origin of the rotated coordinate system [X* — 0,1" — 0 , Z' = 0] at heliographic position (/i,&i) at an 
altitude hp 00 t above the solar surface, the transformation into heliographic coordinates is then given by 

/ l\ + arctan [X f /(Z f 4- hp oot -f R®)] \ 

| bi 4- arctan [Y f /(Z f -f hp 00t + Rq)] J * (A6) 

V v/[A ' 2 + y /2 + W+ h F oot + Rq) t ]/Rq ) 

A semi-circular loop (with radius Rq) can then be defined in this loop plane (A f Y, Z) simply by 

/ X \ ( Rq cos(<^?) \ 

y U o , (A7) 

\ Z / \ Ro sin(<^) / 

where the loop length s — Rq^p can be parametrized with one single parameter {ip). 



B. (APPENDIX) Line-of-Sight Correction on Loop Column Depth 


Column depth of loops with constant cross-section : In order to convert observed emission measures 
EM(x, y) = f n 2 e {x,y,z)dz into local electron densities n e (x } y, z) we need information on the column depth 
f dz . An approximation that is often useful are coronal loops with a constant cross-section w, which can 
be measured from the FWHM as it appears perpendicular to the line-of-sight in the plane of the sky. For 
3-dimensional models of loops parametrized by coordinates ( Xi , j/i, z,*), the angle 4> between the line-of-sight 
and a loop segment can then directly be derived by the ratio of the projected to the effective length of a 
loop segment [i, i+ 1]. 


COS ( \p [iC j , 2/i , 


\/(xi+ 1 - Xj ) 2 + (y<+ 1 - yij 2 
\/{xi+i - x t ) 2 4 {Ui+i - Vi) 2 + (-1+1 - z i) 2 


(Bl) 


yielding the column depth w z along the line-of-sight axis z, 

w 

cos(il>[xi,yi,Zi}) 


w 2 [xi, y,, Ji] - 


(B2) 
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Table 1. Physical Parameters of 30 Coronal Loops (1996- Aug-29, EIT 171 A) 


Loop 

Heliographic 

Azimuth 

Inclination 

Emission 

Loop 

Electron 

Scale 

Scale height 

Filter ratio 

# 

coordinates 

angle 

angle 

measure 

width 

density 

height 

temperature 

temperature 


luh 

£> 

i? 

EM max 

w 

Tie 

A 


j' EIT 


[deg] 

[deg] 

[deg] 

log[cm“ 3 ] 

[Mm] 

[10 9 cm -3 ] 

[Mm] 

[MK] 

[MK] 

1 

247.9, -15.4 

15± 1 

-42± 1 

27.71 

6.8 

2.5 

49 

1.084: 0.07 

1.25± 0.09 

2 

247.7, -15.5 

13± 1 

-49± 5 

27.93 

6.1 

2.3 

50 

1 .09± 0.12 

1.27± 0.10 

3 

247.4, -15.6 

12± 1 

-34± 1 

27.79 

7.4 

2.1 

57 

1.25± 0.13 

1.22± 0.10 

4 

247.3, -15.0 

8± 2 

-26± 1 

27.60 

6.9 

2.6 

57 

1.264: 0.14 

1.27± 0.11 

5 

247.1, -14.9 

8± 3 

-49± 1 

27.34 

6.3 

2.5 

56 

1.22± 0.33 

1.23± 0.10 

6 

247.0, -13.9 

2± 2 

-56± 1 

29.68 

7.0 

3.7 

42 

0.934: 0.08 

1.30± 0.12 

7 

246.1,-14.7 

4± 2 

-36± 1 

27.88 

7.1 

1.7 

62 

1.36± 0.17 

1.1 9± 0.09 

8 

245.2,-14.7 

9± 2 

-33± 2 

27.45 

7.1 

2.2 

57 

1.264: 0.10 

1.23dh 0.06 

9 

247.0,-12.2 

-3± 3 

- 1 2 dh 1 

27.52 

8.1 

2.1 

60 

1.32± 0.53 

1.154: 0.03 

10 

245.8, -13.3 

1± 2 

-23± 1 

27.49 

7.9 

1.4 

53 

1.16± 0.30 

1.214: 0.11 

11 

244.9,-12.2 

1± 3 

-314: 2 

27.11 

6.8 

1.6 

66 

1 .44 ± 0.18 

1.20± 0.09 

12 

247.3,-10.7 

-6± 2 

10± 1 

27.56 

7.8 

2.0 

60 

1 .31 ± 0.07 

1.10± 0.06 

13 

250.6, -10.1 

-9± 7 

114: 1 

27.41 

6.4 

1.8 

60 

1.31 ± 0.17 

1.21± 0.11 

14 

247.1,-9.8 

-8± 3 

10± 1 

27.18 

6.4 

1.7 

33 

0.734: 0.08 

1 .124: 0.05 

15 

247.5, -9.0 

-8db 1 

124: 1 

27.48 

6.8 

2.0 

47 

1.04± 0.08 

1.18± 0.09 

16 

248.1, -8.5 

-11± 3 

25± 1 

27.28 

6.7 

1.4 

53 

1.15± 0.32 

1.24± 0.12 

17 

248.5, -8.0 

-14± 4 

32± 1 

27.04 

7.4 

1.7 

47 

1.03± 0.30 

1.18± 0.04 

18 

249.0, -7.4 

-9± 4 

40± 3 

26.90 

5.6 

1.1 

93 

2.03± 0.41 

1.30± 0.06 

19 

249.6, -6.9 

-18± 3 

39± 1 

27.28 

7.9 

1.3 

65 

1.43± 0.17 

1.18± 0.04 

20 

250.7, -6.6 

-24± 6 

52± 6 

27.43 

7.7 

1.8 

56 

1.22± 0.05 

1 .224: 0.04 

21 

251.5, -6.1 

-21± 6 

584: 3 

27.51 

7.6 

2.1 

45 

0.98db 0.11 

1.15± 0.04 

22 

251.9, -5.7 

-25± 5 

58± 3 

27.32 

7.2 

1.0 

58 

1.27± 0.82 

1.224: 0.12 

23 

259.1, -10.7 

175± 4 

0± 1 

27.61 

8.0 

2.3 

60 

1 .3 1 db 0.28 

1.234: 0.15 

24 

258.6, -11.4 

178± 4 

14± 6 

28.29 

8.7 

1.3 

64 

1 ,41d: 0.22 

1.20± 0.08 

25 

263.2, -16.6 

1784: 5 

43± 1 

27.38 

8.1 

1.8 

42 

0.934: 0.04 

1.18± 0.11 

26 

259.9, -13.8 

1714: 4 

13± 2 

26.95 

5.1 

1.4 

59 

1.304: 0.36 

1.154: 0.08 

27 

259.1, -13.4 

1 77 ± 3 

50± 2 

27.20 

7.1 

1.5 

58 

1.26± 0.30 

1.214: 0.07 

28 

258.0, -13.4 

1794: 2 

274: 2 

27.79 

7.5 

2.6 

57 

1.244: 0.05 

1.124: 0.04 

29 

257.9, -13.7 

185± 2 

36± 1 

27.76 

6.0 

2.4 

60 

1.314: 0.06 

1 .364: 0.07 

30 

257.5, -13.8 

1774: 3 

694: 1 

29.48 

7.4 

1.7 

38 

0.844: 0.10 

1 .164: 0.08 

average 

251.0, -11.8 

-34:10 

7 ±37 

27.61± 0.61 

7.1±0.8 

1.94:0.6 

55± 10 

1.22± 0.23 

1.214: 0.06 
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Fig. 1 — SoHO/EIT Fe IX/X image of active region AR 7986, recorded on 1996-Aug-30, 0020:14 UT, at a 
wavelength of 171 A, sensitive in the temperature range of T e = 1.0— 1.5 MK (top). The greyscales of the 
image is scaled logarithmically in flux, the contours correspond to increments of 100 DN (data numbers). 
The heliographic grid has a spacing of 5°. The filtered image (lower panel) was created by subtracting a 
smoothed image (using a boxcar of 3x3 pixels) from the original image, in order to enhance the loop fine 
structure. 


Fig. 2. — Definition of loop parameters: loop point positions (xi y yi),i = starting at the primary 

footpoint at height h\ = hj 00 t, the azimuth angle a between the loop footpoint baseline and heliographic 
East- West direction, and the inclination angle d between the loop plane and the vertical to the solar surface. 


Fig. 3.— The principle of Dynamic Stereoscopy is illustrated here with an example of two adjacent loops, 
where a thicker loop is bright at time t \ , while a thinner loop is brighest at time ^2- From the loop positions 
(xi, yi) measured at an intermediate reference timet, i.e. t\ < t < ^2 (middle panel in middle row), projections 
are calculated for the previous and following days for different inclination angles d of the loop plane (left 
and right panel in middle row). By extracting stripes parallel to the calculated projections d — 10°,20°,30 3 
(panels in bottom part) it can be seen that both loops appear only coaligned with the stripe axis for the 
correct projection angle i) = 20°, regardless of the footpoint displacement Ax between the two loops. The 
coalignment criterion can therefore be used to constrain the correct inclination angle x ) , even for dynamically 
changing loops. 


pig 4 — Projections of 30 stereoscopically reconstructed loop segments (numbered white curves) are shown 
overlaid on the SoHO/EIT 171 A images (top panels) and the filtered images (bottom panels) of 1998-Aug-29, 
30, and 31. The 30 loop segments were traced from the filter image of Aug 30 (bottom middle panel), while 
the projections on the previous and following day were calculated from the inclination angles t? obtained 
from the dynamic stereoscopy method. Note that the overall magnetic field structure is almost invariant 
during the 3 days, but dynamic changes of the loops produce slight displacements between the calculated 
projections forward and backward in time and the actually observed fine structure. 


Fig. 5. — Three different projections of the stereoscopically reconstructed 30 loops of AR 7986 are shown. 
The loop segments that were traced from the 96-Aug-30 171 A image are marked with thick solid lines, while 
the extrapolated segments (marked with thin solid lines) represent circular geometries extrapolated from 
the traced segments. The three views are: (1) as observed from Earth with /o,&o (bottom right panel), (2) 
rotated to North by &' 0 = 6 0 — 100° (top right panel), and (3) rotated to East by l f 0 = / 0 +97.2° (corresponding 
to -7.2 days of solar rotation; bottom left panel). An EIT 171 A image observed at the same time (-7.2 
days earlier) is shown for comparison (top left panel), illustrating a similar range of inclination angles and 
loop heights as found from stereoscopic correlations a week later. The heliographic grid has a spacing of 5° 
degrees or 60 Mm. 
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Fig. 6. — The positions of the curved coordinate grids of the 30 analyzed loop segements are shown in the 
top right panel, which has the same orientation as the 1996-Aug-30 map shown in Fig. 4 (middle). The 
coordinate grid of loop ^1 is represented with 1-pixel resolution, while only the outer borders and central 
axes are indicated for the other loop segments. - The vertically oriented panels represent the coordinate grids 
of the analyzed 30 loops, stretched out along the loop axis. The top of the panels corresponds to the primary 
footpoint (see positions /i,6i in Table 1). In each panel we show the EIT 171 A flux loop cross-sections 
measured perpendicularly to the loop axes. Successive cross-sections are separated by a distance of 1 pixel 
along the loop axis. The flux associated with each loop is marked with a grey area, obtained by background 
subtraction with a cubic spline interpolation between both sides of the loop cross-section. 


Fig. 7. — The loop-associated EUV flux (F/ oop (s), thick line), the total EUV flux (F to t a i{s) , thin line), and 
the background (F to tai(s) - Fji ux (s ), dotted line) measured along the loop is shown for Loop #1 (left panel). 
A histogram of the relative fractions Fi oop / Ft 0 t a l (integrated over the traced loop lengths) is shown from all 
30 analyzed loops in right panel. 


Fig. 8. — The variation of the loop thickness is shown for the 10 loops with the least confusion by adjacent 
loops (see cross-sections in Fig. 6) as function of the loop length s (left panels). A linear regression fit is 
indicated (solid line in left panels). The average (equivalent) width w is histogrammed for all 30 analyzed 
loops (top right). A loop divergence factor is calculated from the ratio of the average width in the upper 
half and lower half (traced) loop segments (right bottom). Note that most of the loops show no significant 
loop thickness variation. 


Fig. 9.— Left panel: The effect of the variable column depth ic 2 (s) measured parallel to the line-of-sight, 2 
is illustrated as function of the loop length parameter s, for a loop with a constant diameter w. Right panel: 
The effect of the inclination angle d of the loop plane on the inferred density scale height A(tf) is shown. 
Both effects have to be accounted for the determination of the electron density n e (s) along the loop. 


Fig. 10. — The electron density n e (s) (left panels) and the EIT filter-ratio electron temperature T e {s) (right 
panels) as function of the loop length s for the same 10 loops selected in Fig. 8. An exponential model n e (h) 
is fitted, yieling the density scale height A and scale height temperature T x indicated in the left panels. 
The average EIT filter-ratio temperature (obtained from 171/195 A images) is indicated in the right 

panels. A temperature gradient dT e /ds is listed if the gradient is significant. 


Fig. 11.— Statistics of scale height temperatures T x (left top panel), EIT filter-ratio temperatures Tf IT 
(left bottom panel), scatter plot of these two temperatures (right top panel), and EIT temperature gradients 
dT/ds (right bottom panel) for all analyzed 30 loops. 


Fig. 12.- 
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Figure 1 : AR 7986 (Aschwanden et al. 1998) 
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Principle of Dynamic Stereoscopy 
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Figure 3 : AR 7986 (Aschwanden et al. 1998) 
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Figure 7 : AR 7986 (Aschwanden et al. 1998) 
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Figure 10 ; AR 7906 (Aschwanden et al. 1998) 
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